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Abstract 

We study the equilibrium properties of a quantum dot connected to a bulk 
lead by a single-mode quantum point contact. The ground state energy and 
other thermodynamic characteristics of the grain show periodic dependence on 
the gate voltage (Coulomb blockade). We consider the case of almost perfect 
transmission, and show that the oscillations exist as long as the transmission 
coefficient of the contact is less than unity. Near the points where the dot 
charge is half-integer the thermodynamic characteristics show a non-analytic 
behavior identical to that of the two-channel spin-^ Kondo model. In par- 
ticular, at any transmission coefficient the capacitance measured between the 
gate and the lead shows periodic logarithmic singularities as a function of the 
gate voltage. 

PACS numbers: 73.20.Dx, 73.40.Gk 
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I. INTRODUCTION 



The phenomenon of Coulomb blockade of tunneling has recently attracted a lot of inter- 
est, both theoretical and experimental It can be observed, e.g., by measuring conduc- 
tance of a system of two macroscopic leads connected to a small metallic grain by tunnel 
junctions. At low temperature tunneling of an electron into the grain leads to an increase 
of the electrostatic energy of the system by finite amount Ec = e^/2Co, where Cq is the 
grain capacitance. Thus the tunneling conductance becomes exponentially small, with the 
activation energy Eq. One can then add a gate electrode in order to control the electrostatic 
energy 

Here Q is the grain charge; parameter N is proportional to the gate voltage Vg. The 
activation energy is now a function of V^. At the values of the gate voltage corresponding to 
N = n + ^ the energies of states with charges en and e(n + 1) are equal, and the activation 
energy vanishes. Therefore one observes periodic peaks of conductance as a function of the 
gate voltage. 

Recently the Coulomb blockade was observed in semiconductor heterostructures [0. Un- 
like in metallic systems, in a semiconductor device it is often possible to control the barrier 
height by adjusting voltage on additional gate electrodes. In such experiments one can 
study the evolution of the Coulomb blockade as the transmission coefficient T of the tunnel 
barrier changes from to 1. The experiments [0,^ indicate that the increase of the trans- 
mission coefficient leads to the suppression of Coulomb blockade. At T ~ 1 instead of well 
separated peaks weak periodic oscillations of conductance G{Vg) are observed. Experiment 
indicates that the Coulomb blockade disappears at T = 1. On the other hand, in the 
experiment the Coulomb blockade oscillations were observed even when the conductance 
of the junction exceeded e'^/irh. To resolve this contradiction one needs a theory of the 
Coulomb blockade in the regime of strong tunneling, T — 1. 

In this paper we study the Coulomb blockade in a quantum dot connected by a con- 
trollable tunnel junction to a single electrode. Fig. |l|a. Conductance measurements in such 
a system are not possible. However, the Coulomb blockade shows up in the oscillations of 
the equilibrium characteristics of the system, e.g., its ground state energy E, or the average 
charge {Q) of the dot. Experimentally the capacitance C = d'^E/dV^ between the gate and 
the lead can be measured 0. 

Most of the theoretical work on Coulomb blockade is devoted to the case of weak tun- 
neling, when the transmission coefficient of the tunnel barrier is small: T ^ 1. At zero 
temperature, in the limit of very high barrier the charge of the dot Q{N) is quantized in 
units of the elementary charge e, except for the degeneracy points = n + |, where Q 
changes from ne to (n + l)e (solid line in Fig. However, if the small probability of tun- 
neling through the barrier is taken into account, the charge of the dot is no longer a good 
quantum number. As a result the plateaus in {Q{N)) are not horizontal, with the slope 
proportional to T, Ref. 0. This phenomenon is due to the quantum fluctuations of the 
dot charge caused by virtual processes of electron tunneling between the grain and the lead. 
Furthermore, the quantum fluctuations were shown [0 to smear the steps of the average 
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grain charge Q at half-integer values of A^, making {Q{N)) a continuous function (dashed 
line in Fig. |]). 

In this paper we present a theory of the Coulomb blockade near the strong tunneling 
limit T = 1. In Sec. |I|we show that the Coulomb blockade in the system shown in Fig. |l|a is 
described by a one-dimensional (ID) model. This allows us to use the bosonization approach, 
and treat the Coulomb interaction exactly. In agreement with experiment we find no 
contributions which are periodic in Vg in any measurable characteristic of the system at T = 
1. The backscattering on the barrier at T < 1 can be treated in the bosonization approach as 
a small perturbation. In Sec. Q we calculate the first non- vanishing correction to the ground 
state energy of the system E{N), average grain charge {Q{N)), and capacitance C{N). 
These corrections are periodic in A^, with the period corresponding to the change of the 
grain charge by e. In the case of electrons with spin the corrections diverge logarithmically 
at low energies, indicating that the higher-order calculation is necessary. Such a calculation 
performed in Sec. |V| removes the singularities at all values of N except half-integer ones. 
In Sec. we discuss the non-analytic behavior of the thermodynamic characteristics of the 
system a.t N = n + ^ using the analogy [0 between the Coulomb blockade problem and 
the two-channel Kondo model. We argue that the periodic logarithmic singularities in the 
capacitance measured between the gate and the lead should be observed at any value of the 
transmission coefficient T < 1. 

II. ONE-DIMENSIONAL MODEL 

The system we study is shown in Fig. |l|a. The dot is connected to the lead by a narrow 
constriction formed by applying voltage Va to the auxiliary gates. We assume that the width 
of the constricton in its center allows only a single transverse state below the Fermi level. 
In this sense the electron gas inside the constriction is one-dimensional. As the electron 
moves away from the center of the constriction, the the channel becomes wider. Fig. |l|b, and 
the number of transverse modes grows. Since the constriction is formed electrostatically, its 
boundaries are smooth and do not scatter the electrons Thus the constriction creates 
an ideal quantum point contact between two two-dimensional (2D) regions: the dot and the 
lead. 

In the following we will neglect the fact that the dot size is finite, i.e., the infinite 
system shown in Fig. will be considered. The difference between the two systems is 
that an electron entering a finite dot will eventually return back to the lead through the 
constriction. The time r between these two events is determined by the inverse width of the 
discrete energy levels in the grain. In the case of an ideal single-mode junction the width 
is equal to the level spacing e in the dot, and t h/e. An important difference between 
non- interacting and interacting systems is that the latter has another energy scale, Ec- In 
particular, the typical frequency of charge fluctuations is Ec/fi (see, e.g.. Sec. fT|). Since in 
a 2D dot Eq 3> £, the characteristic time Ti/Ec at which the Coulomb blockade develops is 
much shorter than r. In the following we will consider the limit r — *• cxd corresponding to 
the inflnite system. 

An important property of the inflnite system shown in Fig. |I|b is that it is essentially 
one-dimensional. To see that, let us consider wave function '^{x,y) of an arbitrary state 
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penetrating the constriction. Near the center of the constriction the wave function is one- 
dimensionah 



^>,{x,y) = My)e''\ (2) 

where 0o(l/) is the wave function of the ground state for the transverse motion, wavevector 
k is determined by the energy corresponding to y). Since we consider an ideal contact 
characterized by the quantized value of conductance G = e'^/nh, the wave function (|]) does 
not have a scattered component (f)Q{y)e~^^^ . Outside the constriction the wave functions 
have a much more complicated form. In particular, they may be strongly affected by the 
disorder present in the 2D leads. Nevertheless, one can label any wave function by a single 
parameter k. Thus the Hamiltonian of the system of electrons penetrating the constriction 
can be written as 

Ho = J iEo + ^— j ala^dk, (3) 

where aj, is the creation operator for the electron in state \E'fc(x, y), and Eq is the energy of 
the transverse motion corresponding to the wave function (poiv)- 

It is worth noting that Hamiltonian (|^) does not describe the whole system of 2D elec- 
trons. For instance, the electrons with energies below Eq do not penetrate the constriction 
and are not included in Hq. However, since we are interested in electron transport through 
the constriction, the existence of electron states confined in one of the electrodes, and there- 
fore omitted in @, does not affect our results. 

The Hamiltonian (j^) has a one- dimensional form. However, unlike in a usual ID system, 
the density of states is determined by the 2D leads and is therefore energy-independent. 
As a result the physical properties of the system at energy scales of order of Fermi energy 
Ep cannot be described by the ID model. On the other hand, the low-energy properties of 
the system, such as conductance at low voltage and temperature eV,T <C Ep, or Coulomb 
blockade that develops at energy scale Eq -C Ep, can be described by the ID model. 

To consider the Coulomb blockade at T < 1 we have to add scattering potential V{x, y). 
We assume that this potential is localized inside the constriction, where one can use the 
simple ID form of the wave functions. Then the Hamiltonian describing such scattering 
also takes a ID form: 

H' = ^llvik- k')ala,,,dkdk'. (4) 
Here the ID scattering matrix element V{q) is determined as 

Viq) =11 V{x,y)\My)\\-''^^dxdy. (5) 



To complete our ID formulation of the Coulomb blockade problem, we must show that 
the interaction Hamiltonian also has a ID form. Assuming good screening within the 2D 
dot, we will describe the Coulomb interaction by the charging energy (0), with Q being the 
charge inside the dot. In the absence of a tunnel barrier the boundary of the dot is not well 
defined; we will assume it to be at the center of the constriction To find explicit form 
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of the charge operator, we note that there is an obvious relation between Q and the current 
operator: 

g = -^[g,i/o] = J(o), (6) 

where J(0) is the operator of current at point x = 0. The current operator can be obtained 
by integrating the standard expression for the current density over the transverse coordinate 
y. The expression for the current density at a; = is local, and we can use the ID form 
of the wave functions inside the constriction. Then the current operator takes the form 

J = jj{k + k')a\a,,,dkdk\ (7) 

where m is the electron mass. Both Hq and J(0) have ID forms in terms of operators Ofc. 
Hence the charge operator found form Eq. @ is also essentially one-dimensional: 

Q = -i^ [[ d^dkdk'. (8) 
^ 2nJJ k-k' ^ ^ 

Equations (0), (^, (|I|), and (|p present a complete ID Hamiltonian of the system in k- 
representation. It is more convenient to treat this Hamiltonian in coordinate representation, 
which can be obtained by making Fourier transformation to the new ID fermion operators: 



^Ij{x) = -= I akc'^'^dk. (9) 



Unlike the initial 2D wave functions these new operators are completely one- 

dimensional. 

As we already mentioned, we are interested in low-energy properties of the system. 
Thus we can linearize the spectrum of electrons in Eq. (|^) near the two Fermi points, 
and write the fermion creation operators in terms of left- and right-moving fermions: 
iIj{x) = ipiix) + iPr{x). The three parts of the Hamiltonian then transform to 

HQ = hvF j ^l{x){iV - kF)^IJLi^) 

+ kF)Mx)]dx, (10) 

H' = [ V{x)tlj\x)ilj{x)dx, (11) 



_ (Q - eNf 

- 2Co ■ ^^^^ 
Here vf is the Fermi velocity; the ID scattering potential V{x) is obtained from the real 2D 
potential V{x,y) by averaging over the electron density |0o(l/)P in the transverse direction. 
The charge operator (§) takes the form 

Q = \j + i^R{x)i'R{x)] sgnxdx. (13) 

As expected, the operator (|T3|) has the simple meaning of the charge transferred from the 
region x < to the region x > 0. 

To summarize, we established that the Coulomb blockade in a dot connected to a bulk 
electrode by means of a quantum point contact can be treated as a ID problem |TD|. This 
will greatly simplify the following discussion, since we can now bosonize the Hamiltonian 
(p!0|)-([T3|) and treat the quartic in fermion operators interaction term ([T^) exactly. 
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III. PERTURBATION THEORY IN REFLECTION AMPLITUDE 



A. Bosonized Hamiltonian 



The bosonization technique |jTT| is apphcable whenever the system behavior is determined 



by the low-energy properties of the ID electron system. As we already mentioned, the 
typical energy scale at which the Coulomb blockade becomes important is Ec- This energy 
is much lower than the Fermi energy, and the condition of applicability of the bosonization 
approach is satisfied. At low energies the electron system can be viewed as an elastic 
medium. Therefore the bosonized Hamiltonian can be written in terms of two variables: the 
displacement of the elastic medium u{x) and its momentum density p{x). The Hamiltonian 
(p^) of non-interacting electrons takes the form: 

H. = l{^^ + '^m„,4lVui.)r]d.. (14) 

Here m is the electron mass, Uq = mvp/T^h is the electron density. The two fields, displace- 
ment u{x) and momentum density p{x), satisfy the standard commutation relation 

[u{x),p{y)] = ih5{x-y). (15) 

To bosonize the interaction Hamiltonian ([12|) we should find the expression for the charge 
transferred through the point x = 0. This can be done by substitution of the bosonized 
formula for the electron density + '4'r'4'r ~^ —noVu into Eq. (|T^) . Alternatively one 

can just note the obvious relation Q = eno'u(O) between charge and displacement. In either 
case we get 

Hc = Ec[nou{0)-N]^. (16) 

It is important to emphasize that the interaction term (0) is quadratic in bosonic variables, 
and the Hamiltonian Hq + Hq of the system without scattering potential can be diagonalized 
exactly. 

We are primarily interested in the periodic dependence of the ground state energy on 
the gate voltage, i.e., on N. One can easily see that after the transformation 

u{x) u{x) + N/no (17) 

the Hamiltonian Hq + He does not depend on A^. Hence we expect the Coulomb blockade 
to be completely suppressed by charge fluctuations in the absence of scattering potential, 
i.e., at T = 1 (thin line in Fig. 0). To find the Coulomb blockade oscillations at T < 1, one 



has to consider the effect of scattering on the potential of the barrier |jT2|. 

In the bosonization approach only the low-energy properties of the system are considered. 
The scattering on a localized potential is therefore characterized by two constants: the 
amplitudes of forward and backward scattering V{0) and V"(2/c^) given by Eq. (^. In boson 
representation the scattering potential has the form: 

H' = -1/(0)^0 Vm(0) - Y^^hllo cos[27rnoM(0)]. (18) 

nnvF 
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where D is the high energy cut-off (bandwidth). The first term in Eq. (|T8]) does not change 
under the transformation ( |T7|) and does not lead to a dependence of the ground state energy 
on A^. On the contrary, the backscattering term becomes periodically dependent on iV, and 
therefore leads to periodic dependence of all the thermodynamic properties on gate voltage. 
Below we calculate the first non-vanishing periodic correction to the ground state energy 
due to the backscattering. 



B. Perturbation theory for spinless electrons 

To find the first-order correction to the ground state energy one has to calculate the 
average of the cosine term in Eq. (|T8|) over the ground state of the quadratic Hamiltonian 
(|I^), (0). This can be performed along the same lines as in the Debye- Waller theory: 

= --|r|Dcos[27™o(u(0))]e-2"'"o<<"'(°)». (19) 

Here we introduced the reflection amplitude r. (In the first order in backscattering potential 
r = V{2kp)/ihvF.) From Eq. ([T6|) it is obvious that (m(0)) = N/uq. In a ID elastic 
medium the average quantum fluctuation of the displacement ((m^)) = (m^) — (m)^ diverges 
logarithmically due to the low-energy phonons. In our case the phonons with energies below 
Ec are pinned down by the interaction term ([TBI) , see Ref. Therefore the fluctuation of 
u is large, but finite. To find {{u"^)) we will explicitly diagonalize the Hamiltonian Hq + He- 
This is achieved by the transformation 

u{x) = 1 — ^ / UkCos{k\x\ — 5k)dk^ (20) 

rin J TX Jo 



no 

roc 



p{x) = I pkCos{k\x\ — 6k)dk. (21) 



'TT Jo 

Here the phase shift 6k is defined as 

Ec 



arctan 



TrhvpkJ 



) ■ (22) 



In Eqs. ( pOD and ( pT) ) we neglected the odd modes proportional to sin kx because they do 
not contribute to u{0) and therefore are decoupled from both interaction He and scattering 
H'. The new fields Uk and pk satisfy standard commutation relations [uk,Pk'] = ifi6{k — k'). 
In terms of these new fields the Hamiltonian Hq + He takes the form 

Ho + Hc= r (-^ + lmnoivFk)\l] dk. (23) 
JO \2mnQ 2 J 

It follows from Eq. (|20| ) that the contribution of the low-frequency modes to the displacement 
u{0) is suppressed by the factor cos Sk ~ hvpk/Ec. This gives rise to the low-frequency cut- 
off in the logarithmic integral for ((u^(0))). A simple calculation then gives 

^'(0))) = 7r^lnf^V (24) 



27r^nl We, 



c . 
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Here 7 = , with C ^ 0.5772 being the Euler's constant. We can now find the periodic 
correction to the ground state energy using Eq. ([T9|) , 

= -^|r|Eccos27riV. (25) 

As expected the amphtude of oscillations of the ground state energy becomes of the order 
of Eq at |r| ~ 1. 

The average charge (Q) in the dot can now be found using (P^, 

(Q) = eiV - 77^^^ = eiV - ^e|r| sin27riV. (26) 

Civ TT 

The period of oscillations corresponds to the change of the average number of particles in 
the dot by one. At weak refiection, |r| -C 1, the amplitude of the oscillations of charge is 
small (dash-dotted line in Fig. ^). 

Finally, the periodic correction to the capacitance measured between the gate and the 
lead can be found as bC\ = d'^6Ei/dVg. It also exhibits periodic oscillations as a function 
of the gate voltage. 

In the above discussion we completely ignored the spins of electrons. The spinless case 
can probably be realized in an experiment in a high magnetic field. However in the absence 
of magnetic field one should take the spins into account. We will now demonstrate that the 
spin degree of freedom affects the above results dramatically. 



C. Perturbation theory for electrons with spin 



To take into account electron spins we consider a model with two channels corresponding 
to two spin directions: 



Ha 



' [ 2mnQ 2 
He = Ec{nQ[ui{Q) + U2 



-N}\ 



>dx. 



H' 



1 



7r 



'|D{cos[27™oMi(0)] + cos [277^0^2(0)]}. 



(27) 

(28) 
(29) 



The expressions for Hq and H' are obtained by a straightforward generalization of Eqs. (|T4D 
and ( |T^ ) to the two-channel case. (In Eq. ( pQ] ) we neglected the forward scattering terms.) 
The two channels are coupled through the interaction term (^8]) which depends only on the 
sum Mi(0) + ^2(0) representing the total charge brought into the dot. Thus it is natural to 
transform the Hamiltonian to charge and spin modes Mc,s = {ui ± U2)/\/2 (and similarly for 
momentum densities Pc,s)- In the new variables the Hamiltonian (p7D-(p9|) takes the form: 



^0 

He 
H' 



{pi + _^ mnoVp 
I 2mnQ 2 



E, 



c 



(Vmc) + (Vm, 
2 



>dx. 



- — \r\D cos[V2 TTnoUc{0)] cos[-\/2 7rno'Us(0)]. 



vr 



(30) 

(31) 
(32) 
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One can easily calculate the first-order correction to the ground state energy following 
the discussion of the spinless case. An important difference between these two cases is that 
now we have two modes, and only one of them is pinned down by Coulomb interaction (^TJ). 
Therefore the quantum fluctuation {(u^)) of the displacement in the spin channel diverges 
logarithmically. This leads to a strong suppression of the oscillations of the ground state 
energy. The amplitude of these oscillations is \r\^/Ece, where e is the low energy cut-off 
of the order of the level spacing in the dot. In the limit of large dot {e — > 0) the first- 
order correction to the ground state energy vanishes. Thus in order to obtain the Coulomb 
blockade oscillations one has to perform the calculation up to the second order in barrier 
potential. 

The second-order correction to the ground state energy can be presented in the form 

1 r°° 

6E2 = -Im / {H\t)H'{0)) dt. (33) 
Ti Jo 

From the explicit form (|32| ) of the perturbation H' it follows that the correlator {H'{t)H'{0)) 
factorizes into charge and spin parts. The spin part is easily calculated: 

(cos[v^7™oM.(0,t)]cos[y2 7moM.(0,0)]) = (34) 

The slow decay of the correlator at large t is due to the low-frequency modes. In the charge 
channel the low-frequency components of u{0) are suppressed. As a result at t h/Ec the 
charge part of the correlator {H'{t)H'{0)) saturates, 

(cos [v^7™oMc(0, t)] COs[V27™oMc(0, 0)]) 

= COS ttN ^-^ sm vrA/ . (35 

TiD \ AElt^ I ^ ' 

The substitution of Eqs. (|^) and (|35|) into the expression for the second-order correction 
(PB|) gives the integral which diverges logarithmically at large t. The divergence can be cut 
off at t ~ /i/e (or at t ~ h/kBT if the correction to the free energy at a finite temperature 
is being calculated). The result has the form 

5E2 = -^\r\^Ecln (^^) cos^nN. (36) 

In the limit of large dot Eq/b — > 00, and the second-order result diverges. This indicates 
that the terms of higher orders in |r| should be taken into account. 



IV. HIGHER-ORDER CALCULATION OF THE GROUND STATE ENERGY 

To proceed with the higher-order calculation we will first simplify our Hamiltonian (|30D - 
(p^). Since the logarithmic divergence arises at small energy scales E <C Ec, we do not 
have to treat the charge fluctuations exactly. At such low energies the charge fluctuations 
are suppressed by the interaction term, and one can replace cos[v^ 7rnoUc(0)] in Eq. ( P^ by 
its value averaged over the unperturbed ground state. After this simplification the charge- 
related part of the Hamiltonian completely decouples and can be excluded. 
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Another simplification is made possible by the fact that the barrier potential depends 
only on the spin mode displacement at x = 0. Therefore the odd elastic modes proportional 
to sinkx are not coupled with H' and can be excluded. To this end we will change the 
variables: Ue,o{x) = [usix) ± Us{—x)]/\/2, and similarly for the momentum densities Pe,o(2^)- 
Thus we arrive at a Hamiltonian 



Pe 



H' 



\2mnQ 
8-fEcD 



+ 



{Vu,f\dx, 



1/2 



|r| cos7rA^cos[7™oMe(0)]. 



(37) 
(38) 



The bandwidth D here should be taken of the order of the charging energy Ec-, because 
only at such small energies the charge fluctuations can be neglected. 

The Hamiltonian (|37D, (^) is very similar to the Hamiltonian (14), ([T8| ) of non- 
interacting electrons in the presence of scattering potential. The important difference is 
that the cosine in Eq. (pH]) has twice smaller argument than the one in Eq. (|18]). The latter 



_I8]) has twice smaller argument than the one in Eq. 
represents a product of two fermion operators ■?/'^(0)^/'(0). Similarly, Eq. (|38D can be inter- 
preted as a sum of two fermion operators: '?/'^(0) + V'(O)- To support this observation we will 
re-write the Hamiltonian in terms of a new bosonic field 



-'KnQUe{\x\) + 



sgnx 
TinQ Jo 



{x')dx' . 



(39) 



Unlike the old variables Ue{x) and Pe{,x) defined at a; > 0, the new field $(a;) is defined on the 
whole X-axis and has the same number of degrees of freedom. The commutation relations 
for the new variables have the form: [$(a;), $(y)] = ^vr sgn [x — y). 
In terms of the field $ the Hamiltonian (|37f), (BSl) takes the form: 



H' 



Tivp 
A-n 



[V^{x)fdx, 



S'-fEcD\^/'^ 

— Irl cos7rA^cos$(0). 



(40) 
(41) 



Expression (^0]) obviously coincides with the well-known form of the bosonized Hamil- 
tonian of a ID gas of non- interacting right-moving electrons (see, e.g., Ref. ||ll|). In this 
' D / 2'Khv p e^'^'^^'^ is identified as the electron annihilation operator ip^x). 



case the operator 

Consequently, the Hamiltonian (^01), (^Tj) can be de-bosonized to the following form: 

Hf = r [ikhlh, - \{hi + h,)] dk, (42) 

J —CO '- -' 

where bf^ and bl are the new fermion operators in ^-representation, electron energy = hvpk 



is measured from the Fermi level, and A = y2'yhvpEc/TT^ \r \ cosnN. 

The Hamiltonian (^2]) contains a term linear in fermion operators b/^ and b].. If it is 



treated perturbatively, one obtains the result (pq). However, it is possible to take into 
account all the higher terms of the perturbation theory. This can be done by transforming 
the Hamiltonian (^2]) to a quadratic form. One such transformation based on the Jordan- 



Wigner transformation to a spin chain was suggested by Guinea ||1J]. We will use a simpler 
t r ansf ormation : 
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bk = {c + c^) Ck. (43) 

One can easily check that if c and Ck are fermion operators, the operators bk have correct 
anticommutation relations. After this transformation we get a Hamiltonian quadratic in 
fermion operators 

Hg = j^Ji^kcW - A [4 (c + ct) + (c + c^) cj} dk. (44) 

The Hamiltonian Hq is very similar to the Hamiltonian of a resonant impurity at the 
Fermi level. Unlike the latter, Hq does not conserve the number of particles and should be 
diagonalized by a Bogoliubov transformation. As a result we get the diagonal form 

Hq = E + a {CIC, + ClC,) dk, (45) 

where E is the ground state energy of our Hamiltonian. The two branches of the excitation 
spectrum correspond to some linear combinations of particle and hole states. One of the 
braches is not affected by coupling to the impurity level: Cp, = {c^ + cl^)/v^. The other 
branch has some admixture of operators and c. 



r - ~ ^-^ / ^^^^ (r I A 



where the principal value of the integral is assumed. The parameter F = AnX'^/hvp has the 
meaning of the width of the resonant level. 

The correction SE to the ground state energy of the Hamiltonian Hq can be found, e.g., 
by averaging Eq. (|45|) over the unperturbed ground state, 

SE = - UClC,),dk. (47) 





The resulting integral over k is logarithmically divergent at the upper limit due to the second 
term in the right-hand side of Eq. (^61). However, as we already mentioned, the bandwidth in 
our Hamiltonian should be Ec- Unlike in the case of the perturbation theory, the logarithmic 
integral now has an intrinsic cut-off F at low energies. This low-energy cut-off is due to the 
higher order terms in A. As a result the correction to the ground state energy is now finite, 
5E = — (F/27r) \n{Ec/T). In our original notations this result has the form: 

SE = -^\r\^Ec\n ( -—1——] cos^ vrAT. (48) 

TT-^ \ \r\^ COS^ 7TI\ J 

To summarize, the chain of transformations has lead us from bosonized Hamiltonian 
(P7|)-(PP|) to simple form (|1^). The transformations are exact at low energies. They uncover 
the low-energy cut-off F for the logarithmic divergence of the second-order perturbation 
theory. 
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It is important to note that F vanishes at half-integer values of A^. At these points the 
logarithm in Eq. (^SD diverges. Due to the pre-factor, correction 5E is not divergent, but still 
has a non-analytic behavior at N = n + ^. The non-analiticity shows up in the capacitance 
C = d'^E/dVg measured between the gate and the lead: 

SC{N) = ^\r\'' f3^ Ec ln( . \ ^\ os27riV. (49) 
TT V rl cos"' vrA' J 

Here (3 is the parameter controlling the relation between and the gate voltage, A^ = /3Vg] 
its value is determined by the system geometry. One can easily see that the capacitance is 
logarithmically divergent at half-integer A^. The nature of these singularities is discussed in 
the next Section. 

V. ANALOGY TO THE TWO-CHANNEL KONDO MODEL 

The logarithmic divergence of the capacitance ( ^U] ) at half-integer A^ was found for the 
case of strong tunneling. It is instructive to compare it with the similar divergence in 
capacitance in the weak tunneling case: 

*CW = -2^^Ec^exp(j^)l„(^), (50) 

Here t is the transmission amplitude; expression (|50|) is written for the vicinity of the point 
N = ^. Both expressions (^) and (0) predict a logarithmic singularity at A^ = |, with 
the factors in front of the logarithm being of the same order at |t| ~ |r| ~ 1. It is therefore 
natural to conjecture that the logarithmic divergencies of capacitance exist not only in the 
limiting cases of weak and strong tunneling, but at any value of the transmission coefficient. 

To support this idea we shall first outline the arguments leading to the divergence (|50D 
of capacitance in the weak tunneling case. The solution ^ was based on the mapping of the 
Coulomb blockade problem onto an anisotropic multichannel Kondo model. At A^ close to 
I one can consider the perturbation theory in tunneling amplitude and neglect all the terms 
involving virtual states with a charge different from and 1. This restriction on the possible 
charge states is due to the large charging energy associated with all other states. Thus all 
the relevant terms of the perturbation theory are constructed in such a way that first an 
electron tunnels through the barrier from left to right, changing the dot charge from to 1, 
then another electron tunnels from right to left returning the dot to the state with Q = 0, 
then one more electron tunnels from left to right leading to Q = 1, and so on. One can note 
that the same structure of the perturbation theory takes place for the Kondo model with 
anisotropic coupling J±{a^S^ + a^S^). Instead of the two types of electrons, left and right, 
we now have two other types, spin-up and spin-down. Furthermore, each electron scattering 
on the impurity fiips its own spin, e.g., from up to down, and the spin of the impurity, down 
to up. This means that the next electron scattered on the impurity has to fiip its spin from 
down to up, then from up to down, etc. This leads to the same structure of the perturbation 
theory as in the Coulomb blockade problem. 

A small deviation of the system from the point N = ^ gives rise to the energy difference 
between the states with Q = and Q = 1. It is completely analogous to the effect of 
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magnetic field h in the Kondo problem. Thus the capacitance of the system C ~ d'^E/dN'^ 
is analogous to the magnetic susceptibility of the Kondo impurity Xi = d'^E/dh'^. The latter 
is inversely proportional to the Kondo temperature, which leads to the exponentially large 
factor in Eq. (|50D . 

Finally, the presence of the real spin of electrons (conserved in the tunneling process) 
should be interpreted as a "color" for the electrons in the Kondo model. Thus the spin-| 
case maps onto the two-channel Kondo model. The latter is known |T^-0 to exhibit some 



non-Fermi-liquid properties. These include the logarithmic divergence of the susceptibility 
at h = 0, resulting in the logarithmic singularity in Eq. (pD[). 

The result ( |50D for the capacitance in the weak coupling regime was obtained using the 
exact solution of the Kondo model |]15|. Another approach based on the renormalization 



group treatment |T9| allows one to find some low-energy properties of the system at any cou- 
pling strength. The main idea of this technique is that at low energies the effective coupling 
constant in the Kondo model grows, and the system approaches the strong coupling fixed 
point. This fixed point is stable, and the low-energy properties of the system are determined 
by the leading irrelevant perturbation. Therefore the low-energy behavior is universal, i.e., 
independent of the initial conditions. The stable fixed point for the multichannel Kondo 
problem was studied in detail by conformal field theory methods []TB[. In particular, the 



logarithmic behavior of the susceptibility in the two-channel Kondo model was rederived. 

To apply this method to the Coulomb blockade problem one should first find the fixed 
points. There are two obvious fixed points: T = and T = 1. As we already discussed, the 
weak tunneling fixed point T = corresponds to the weak coupling fixed point in the Kondo 
model. This fixed point is therefore unstable. It is then natural to assume that the strong 
tunneling fixed point T = 1 maps to the strong coupling fixed point of the Kondo model. 
To test this hypothesis one should first show that the T = 1 fixed point is stable. This 
can be done by calculating the scaling dimension of the perturbation H' in the Hamiltonian 
The correlator {H'{t)H'{0)) can be found from Eqs. (|3|) and At half-integer 



values of (corresponding to zero magnetic field in the Kondo model) it decays as 1/t 
Thus the scaling dimension of the perturbation if' is |, and the fixed point is stable. 

The stability of the T = 1 fixed point suggests that at half-integer A^ the system with any 
value of transmission coefficient T approaches the strong tunneling limit. We therefore can 
argue that periodic logarithmic divergences of the capacitance (^) and ( |50[ ) are a general 
property of the system and should be observed at any value of the transmission coefficient. 
In particular, the divergence of capacitance at weak tunneling ( |50| ) can be interpreted as a 
consequence of the divergences ( ^Of ) in the strong tunneling limit. 

The last statement can be tested by calculation of the so-called Wilson ratio in the 
strong tunneling limit. Indeed, since the weak-tunneling properties of the Coulomb blockade 
problem were found by mapping to the Kondo model, one expects the correction to the 
heat capacity to have a non-analytic temperature dependence 5cj ~ ksT \q.{Ec / ksT) , and 



its ratio to the correction to the susceptibility to be universal iTqJTql . If the low-energy 
behavior of the weak tunneling model is controlled by the strong tunneling fixed point, the 
same must be valid for the model (pOD-(p^). A straightforward calculation of the correction 
to the heat capacity in the second order in H' gives 



|2 



5q = "^^kBTlniEc/kBT). (51) 



c 
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Taking into account the relation h = 2Ec{N— |) between the effective magnetic field h in the 
Kondo model and the deviation of from a half-integer value, we can find the susceptibility 
from Eq. (||) as 6xi = {2Ecy^d'^6E2/dN\ The resuh is 

6x^ = ^^HEc/kBT). (52) 

The ratio ksTSxi/Sci now takes the universal value 2/7r^, in agreement with the theory of 
the two-channel Kondo model . 

Finally, it is instructive to discuss the relation between our quadratic Hamiltonian ( PI) 
describing the low-energy properties of the Coulomb blockade problem and the similar Hamil- 
tonian for the Toulouse limit of the two-channel Kondo model. The latter was obtained by 



Emery and Kivelson |17| and has the following form: 



H=ivf / tp^x) — - — ax 

J —CO Ox 

[^t(o) + ^(o)][dt_rf] + /,(rftrf_i/2). (53) 



v27ra 



Unlike in the Hamiltonian (|53|) , the second part of Eq. (^41) is proportional to the weak 
magnetic field: X oc h. However, one can still bring the Hamiltonian (|53D to the form 
identical to (^4|). A linear transformation of fermion operators analogous to ( ^61) allows 
one to absorb the second term of Eq. (|53D into the kinetic energy term. After such a 
transformation the third term of Eq. (|53| ) takes the form, which is identical to the second 
term of Eq. ( ^^ at low momenta. Therefore at low energies the Hamiltonians (|4^ ) and ( |53D 
are equivalent. This proves our conjecture that the strong tunneling fixed point is identical 
to the strong coupling fixed point of the two-channel Kondo model. 



VI. CONCLUSION 

In this paper we have studied the Coulomb blockade in the limit of the barrier transmis- 
sion coefficient T close to unity. The Coulomb blockade oscillations of the dot charge and 
capacitance persist as long as T < 1. As T approaches unity, the sharp peaks in the system 
capacitance transform into a weak oscillation with periodic logarithmic divergences at the 
points where the dot charge is half-integer. 

The analogy between the Coulomb blockade and Kondo problem discussed in Sec. |V] 
allowed us to conclude that the logarithmic divergences in capacitance should be observed 
at any value of T. Of course, the exact calculation of the capacitance in terms of the 
transmission and reflection amplitudes is possible only in the limits T ^ 1 and 1 — T ^ 1. 
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FIGURES 



FIG. 1. (a) Schematic view of a quantum dot connected to a bulk 2D electrode. The dot is 
formed by applying negative voltage to the gates (shaded). Solid line shows the boundary of the 
2D electron gas (2DEG). Electrostatic conditions in the dot are controlled by the gate voltage 
Vg. Voltage Va applied to the auxiliary gates controls the transmission coefficient T through the 
constriction, (b) Constriction between two 2D regions. Inside the constriction the wave functions 
have ID form (^). 

FIG. 2. The average charge Q of the dot as a function of dimensionless gate voltage at 
different values of transmission coefficient: T = (solid line), T <C 1 (dashed line), 1 — T <C 1 
(dash-dotted line), and T = 1 (thin line). 
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